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ABSTRACT 

By using the Rayleigh- Ritz method of minimizing the total 
potential energy of a structural system, combined load (me- 
chanical or thermal load) buckling equations are established 
for orthotropic rectangular sandwich panels supported under 
four different edge conditions. Two-dimensional buckling in- 
teraction curves and three-dimensional buckling interaction 
surfaces are constructed for high-temperature honeycomb- 
core sandwich panels supported under four different edge con- 
ditions. The interaction surfaces provide easy comparison of 
the panel buckling strengths and the domains of symmetrical 
and antisymmetrical buckling associated with the different 
edge conditions. Thermal buckling curves of the sandwich 
panels also are presented. The thermal buckling conditions 
for the cases with and without thermal moments were found 
to be identical for the small deformation theory. In sandwich 
panels, the effect of transverse shear is quite large, and by ne- 
glecting the transverse shear effect, the buckling loads could 
be overpredicted considerably. Clamping of the edges could 
greatly increase buckling strength more in compression than 
in shear. 

KEYWORDS Sandwich panels; Mechanical buckling; 
Thermal buckling; Buckling interaction surfaces; Buckling 
interaction curves. 

NOMENCLATURE 

A Fourier coefficients of trial function for w, in. 

Aj • extensional stiffnesses of sandwich panel, 

' 3 lb/in. 

a length of sandwich panel, in. 

Qo edge length of square sandwich panel, in. 

a 'j „ coefficients of characteristic equations 

mnfc/ 

R Fourier coefficients of trial function for 7 **, 

in./in. 

b width of sandwich panel, in. 

n Fourier coefficients of trial function for 7 yZ) 

w mn 

in./in. 


Da 

Dqx 1 Dqv 

bending stiffnesses of sandwich panel, in-lb 
transverse shear stiffnesses in xz, yz planes, 
lb/in. 

flexural stiffness parameters, y/D\\D 22 , in-lb 
Young’s moduli of face sheets, lb/in 2 
effective transverse shear moduli of 

D* 

E Xi E y 

Gcxx^Gcy* 

honeycomb core, lb/in 2 

G X y 

shear modulus of face sheets, lb/in 2 

h 

depth of sandwich panel = distance between 


middle surfaces of two face sheets, in. 

h c 

depth of honeycomb core, h c = h — t Sf in. 

I a 

moment of inertia, per unit width, of a face 


sheet taken with respect to horizontal 
centroidal axis of the sandwich panel, 


Is = \t s h 2 + in 4 /in. 

i 

index, 1, 2, 3, ■ • 

3 

index, 1, 2, 3, • • • 

k 

index, 1, 2, 3, • • • 

k x ,ky 

compressive buckling load factors in x- and 
N a 2 AU 2 

y-directions, k x — ** i r 2 D* ’ 


for a = constant 

hxy 

shear buckling factor, k X y — ^2 

0 = constant 

kxi hy 

modified compressive buckling load factors in 

.. . 7 A/jQq L 6 L. 

x- and y-directions, k x = ^Tp+ 


= Jfc v ^, for ab = a 2 0 = constant 

7 iD 

h X y 

modified shear buckling load factor, k X y 

= fc xy |, for ab = al = constant 

l 

index, 1, 2, 3, • •• 

M X , My 

bending moment intensities, (in-lb)/in. 

Mxy 

twisting moment intensity, (in-lb) /in. 



thermal moments, (in-lb)/in. 

77i number of buckle half waves in x-direction 


N x ,N y 

normal stress resultants, lb/in. 

N xv 

shear stress resultant, lb/in. 

n t n t n t 

i¥ x i 1 y * iy xy 

thermal forces, lb/in. 

n 

number of buckle half waves in y-direction 

Qx,Qy 

transverse shear force intensities, Ib/in. 

T 

temperature, °F 

Ter 

critical buckling temperature, °F 

t. 

thickness of sandwich face sheets, in. 

V 

total potential energy of sandwich panel, 
in- lb 

Vl 

strain energy of sandwich panel, in- lb 

V 2 

work done by external forces, in-lb 

u,v,w 

middle surface displacement components in 
x-, y-, and 2 -direction, in. 

x,y,z 

rectangular Cartesian coordinates 

a xi a yi a xy 

coefficients of thermal expansion, in/in-°F 

Ixx i lyz 

transverse shear strains in xz- and 
j/ 2 -plane, in./in. 

c 

numerical coefficient of Nj in a 

V 

numerical factor in buckling equation, and 
associated with an edge condition 

t 

numerical coefficient of Nj in a]^ nkl 

Vxy, Vyx 

Poisson ratios of face sheets, also used for 
those of sandwich panel 


INTRODUCTION 


The structural components of hypersonic flight vehicles 
(e.g., spacecraft, rockets, reentry vehicles, aircraft, etc.) are 
subjected to hyperthermal loadings caused by hostile aerody- 
namic heating during ascent and reentry, or caused by solar 
radiation during spaceflight. These structural components 
have to operate at elevated temperatures and, therefore, are 
called hot structures. Because of nonuniform heating (which 
is magnified by the cooler substructural frames that act as 
heat sinks) and the mechanical structural constraints, severe 
thermal stresses could build up in those hot structures. Ex- 
cess thermal loading may induce material degradation, ther- 
mal creep, thermal yielding, thermal buckling, thermal crack 
fracture after cool down, etc. Any disruption of the surface 
smoothness of these structures (e.g., metallic thermal protec- 
tion system (ref. 1) or hypersonic aircraft engine inlet struc- 
tures (refs. 2, 3), etc.) caused by the above failure modes, 
especially thermal buckling, could disturb the flow field, cre- 
ating hot spots that could cause very serious consequences to 
the structures. Thus, the thermal load is a key factor in the 
design of hot structures. Reference 1 discusses various design 
concepts of both hot and cryogenic structural components 


for hypersonic flight vehicles. The potential candidates of 
high-buckling-strength hot-structural panels (fabricated with 
superalloys) for hypersonic aircraft applications are tubu- 
lar panels, beaded panels, truss-core sandwich panels, hat- 
stiffened panels, honeycomb-core sandwich panels, etc. (refs. 
4, 5). The combined-load buckling behavior of tubular panels 
was studied by Ko et al. (ref. 4) extensively both theoretically 
and experimentally. The compressive buckling characteristics 
of the beaded panels were investigated by Siegel (ref. 5). 

Recently Ko and Jackson (ref. 6) and Percy and Fields 
(ref. 7) studied the compressive buckling behavior of a hat- 
stiffened panel designed for application to the hypersonic 
aircraft fuselage skin panel. Furthermore, Ko and Jack- 
son conducted simple analysis of thermal behavior (thermal 
buckling of face sheet) of a honeycomb- core sandwich panel 
(ref. 8) and compared the relative combined-load buckling 
strengths of truss-core and honeycomb-core sandwich panels 
(ref. 9). They also investigated the effect of fiber orientation 
of a metal-matrix face sheet on the combined-load buckling 
strength of honeycomb- core sandwich panels (refs. 10, 11). 
Most of the past mechanical buckling analyses of sandwich 
panels (refs. 4-7 and 9-12) and flat plates (refs. 13, 14) were 
conducted for simply supported edge conditions because the 
analysis was mathematically less involved. For the case of 
clamped edge conditions, Green and Hearmon (ref. 15) stud- 
ied combined loading stability of plywood plates, and Smith 
(ref. 16) considered only pure shear buckling of the plywood 
plates. Kuenzi, Erickson, and Zahn (ref. 17) considered also 
shear stability of flat panels of sandwich construction. The 
workers cited here ignored the transverse shear effect in their 
analyses. King (ref. 18) analyzed the stability of clamped 
rectangular sandwich plates subjected to in-plane combined 
loadings, taking into account the rotational effect of the sand- 
wich core. A less-compact displacement function (that could 
be reduced to a simpler Green and Hearmon displacement 
function (ref. 15)) was used, resulting in a very complicated 
expression for the potential energy of the sandwich system. 
Most of the past thermal buckling analysis was done on single 
plates (refs. 19-22) or laminated composite plates (refs. 23- 
27), for which the transverse shear effect may be neglected. 
In the actual application of hot structural panels, most panel 
boundary conditions are closer to the clamped edges rather 
than the simply supported edges. Therefore, this paper will 
consider the combined-load mechanical and thermal buckling 
of sandwich panels under different types of edge conditions by 
taking into account the transverse shear effect, and will com- 
pare the buckling interaction curves and surfaces for different 
edge conditions. 

DESCRIPTION OF PROBLEM 

Figure 1 shows the geometry of a rectangular honeycomb- 
core sandwich panel having identical face sheets. The exten- 
sional and bending stiffnesses of the panel will be provided by 
the two face sheets only, and the transverse shear stiffnesses 
by the honeycomb core only. 


2 




Fig. 1. A honeycomb- core sandwich panel. 

This type of sandwich panel, when fabricated with high- 
temperature alloy (e.g., titanium alloy), becomes the so- 
called hot structure and could be a potential candidate for 
hypersonic aircraft structural applications (ref. 1). Figure 2 
shows the sandwich panel subjected to combined compres- 
sive and shear loadings in its middle plane. The conventional 
Rayleigh-Ritz method of minimizing the panel total poten- 
tial energy will be used in the combined-load buckling anal- 
ysis, accounting for the transverse shear effect (fig. 3). The 
sandwich panel wiU be supported under four different edge 
conditions: 

• Case 1: Four edges simply supported (4S edge condition) 

• Case 2: Pour edges clamped (4C edge condition) 

• Case 3: Two sides clamped, two ends simply supported 
(2C2S edge condition) 

• Case 4: Two sides simply supported, two ends clamped 
(2S2C edge condition) 

where the sides and ends are parallel to the i- and y-axes, 
respectively. 



The problem is to study the effects of the panel edge con- 
dition and the panel aspect ratio on the combined-load buck- 
ling behavior of the sandwich panel. Case 1 has already been 
solved and was published in reference 9. For completeness, 
however, some key equations for Case 1 will be repeated in 
this paper. 

GOVERNING EQUATIONS 
Constitutive Equations 

Fbr the classical orthotropic thick plate theory, the 
thermoelastic constitutive equations for membrane forces, 
moments, and the transverse shear constitutive equation may 
be written as (fig. 2) 



Fig. 2. Forces and moments acting on differential element of 
a sandwich panel. 
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For the sandwich panel whose extensional and bending 
stiffnesses are provided only by the two identical face sheet s, 
and the transverse shear stiffnesses only by the honeycomb 
core, the extensional and the bending stiffnesses A,} 
in equations (1) and (2), and the transverse shear stiffnesses 
{Dqx, Dq v } in equation (3) may be written as 
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where the 2 in front of {t,, /,} in equation (4) is associated 
with two identical face sheets. 

The thermal forces {Nj t Ny t N X y} and the thermal mo- 
ments Af^,} appearing in equations (1) and (2) 

are defined by: 
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where t = 1, 2 are associated, respectively, with the lower 
and the upper face sheets, and [ ] i (i = 1, 2) implies that the 
material properties are associated with temperature (t = 
1, 2). The thermal force and thermal moment contributions 
from the honeycomb core were neglected. 

Energy Equations 

Based on the small deformation theory, the strain energy 
Vi of the heated sandwich panel may be written as (refs. 23, 
24, 26, 27) 
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For the buckling problem, the work done by the in- plane 
orces to produce transverse deflection is given by 

, ( 1 dxdy (8) 


+N U 


The total potential energy V of the sandwich panel is then 
V = Vi + V 2 . 

Panel Boundary Conditions 

The sandwich panel is to be supported at its four edges 
tinder the following four cases of boundary conditions: 

For mechanical buckling: 

Case 1. 4S edge condition: x = 0, a : w = M x = 7 y * = 0; 
y = 0,6 : w = Af y = 7 XX = 0 

Case 2. 4C edge condition: x = 0, a : w = = 7 xr = 

7y* “ 0; y = 0, 6 : u; = ^ = 7x* = 7y* = 0 

Case 3. 2C2S edge condition: x = 0, a : u; = M x — 7 yz = 0; 

y = 0, 6 : w = = 7 XZ = 7 yx = 0 
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Case 4. 2S2C edge condition: x = 0,o : w = — j xz — 

7 yz = 0; j/ = 0, 6 : tx; = M y = = 0 

For thermal buckling: 

In addition to the above boundary conditions, the following 
edge condition is to be imposed: 

x = 0 , o : u = v = 0 ; y = 0 , 6 : u = v = 0 


BUCKLING ANALYSIS 

The conventional Raleigh- Ritz method of minimization of 
total potential energy V will be used in the buckling analysis. 

Panel Deformation Functions 

For an eigenvalue solution via the Rayleigh-Ritz method, 
the trial functions for the sandwich panel deformation 
{w,7**,7y*}» satisfying the boundary conditions, may be ex- 
pressed in the following double Fourier series: 

Case 1. 4S edge condition (ref. 9) 
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Case 2. 4C edge condition (ref. 15) 
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Case 3. 2C2S edge condition (ref. 15) 
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Case 4. 2S2C edge condition (ref. 15) 
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Uniform Temperature (Zero Thermal Moments) 

It is well known that thermal stresses are not caused by ex- 
ternal loads but are the consequences of restrained thermal 
distortion. The intensities of thermal stresses will change 
when a structure is deformed; therefore, the thermal stress 
levels are the functions of strains. In classical thermal buck- 
ling of a structural panel (under uniform temperature rise 
with constrained edges), the first-order lateral deflections of 
the panel will cause only second-order small changes in the in- 
plane strains (thus, thermal stresses) at the onset of thermal 
buckling (ref. 28). In mechanical buckling, however, the ex- 
ternal loads are held constant during buckling. If the second- 
order effect is neglected, then the in-plane thermal loads may 
be considered constant during thermal buckling. Thus ther- 
mal buckling problems would be equivalent to mechanical 
buckling problems; therefore, the conventional methods of 
structural stability analysis may be applied to the thermal 
buckling analysis. 

The buckling equations will be developed first for the 
mechanical budding under the combined loading condition 
{N x -> -N X , Ny - —Ny> N X y -4 -N X y}. The resulting 
mechanical buckling equation could then be applied directly 
to the thermal buckling of the sandwich panels with con- 
straint edges under uniform panel temperature (i.e., {N x = 
-Nj, Ny = —Ny , Nxy = -Njy}, {Mj, Ml, Mj v } = 0, 
u = v = w = 0). 

Characteristic Equations in Terms of Load Factors 


After substituting the trial deformation functions (eqs. (9) 
through (20)) into the energy equations for Vi (eq. (7)) and 
V 2 (eq. (8)) (signs of forcing functions reversed), and after 
performing the double integrations, the components of Vj and 
y 2 may be calculated for different indicial conditions under 
different panel edge conditions (ref. 29). 

Substituting Vi and V 2 expressed in terms of Fourier coef- 
ficients A mn , B mn and C mn (ref. 29) into V, and then min- 
imizing V with respect to each of {A mn , B mn , and C mn ) 
according to the Rayleigh-Ritz principle, 


dV dV dV 
dA mn ~ 8B mn dC mn 


( 21 ) 
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one obtains three homogeneous simultaneous equations (i.e., 
characteristic equations) for each indicial set of {m, n}. 
Those three equations may be combined into one character- 
istic equation containing only the panel deflection coefficient 
Akt . Namely, 


OO 00 A M 
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A kt = 0 


(22) 


where the stiffness/ geometry parameter M mn kt appearing in 
equation ( 22 ) is defined as 
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where the coefficients o^ nki (i,j = 1 , 2, 3) are the partial 
functions of {Aj» A^x, A^y, m » n > a > **» ^y} an< * ^eir 
functional forms vary with indicial and edge conditions. For 
the 4S, 4C, 2C2S, and 2 S 2 C edge conditions there are 1 , 9, 
3 , and 3 different sets of a£ nW , respectively, and these are 
defined in reference 29. 

In equation (23) 77 is a numerical parameter, and in equa- 
tion ( 22 ) 6 m nki is a special delta function which is nonzero 
only when m ± A: = odd, and n ± l = odd. 


The characteristic equation (22) forms a system of infi- 
nite number of simultaneous homogeneous equations (i.e., 
one infinite series equation for each set of {m,n} values). 
Because there is no coupling between the even case (sym- 
metric buckling) and the odd case (antisymmetric buckling) 
(ref. 9 ), those simultaneous homogeneous equations may be 
divided into two groups that are independent of each other: 
one group in which m ± n is even, and the other group in 
which m ± n is odd (refs. 9, 13, 14). 

For the deflection coefficients Akt to have nontrivial so- 
lutions under the assigned values of {k x ,k yi ^}, the deter- 
minant of coefficients of unknown Am of the simultaneous 
homogeneous equations generated from equation ( 22 ) must 
vanish. The largest eigenvalue 1 / kx V thus obtained will give 
the lowest shear buckling load factor k xy for given values 
{Jfc*, ikyi |[}* When the transverse shear effect is neglected 

(eqs. ( 22 ), (23)), {k x ,k y ,k xy } are a function only of §, and 
independent of panel size. However, they will become panel- 
size dependent if the transverse shear effect is considered. 

In the actual eigenvalue computations, the determinants of 
order 12 were found to give sufficiently accurate eigenvalue 
solutions (ref. 9). These determinants are given in reference 
29 for the cases m ± n = even (symmetric buckling) and 
m ± n = odd (antisymmetric buckling) for different edge 
conditions. 

Characteristic Equations in Terms of Temperature 


Case 1. 4S edge condition 
’ 7] =32 

c mnkt 

. 0mn «- (m* - fc 2 )(n 2 " *“) 
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Case 2 . 4C edge condition 
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For thermal buckling, the main objective is to find the 
buckling temperature, T c r , rather than thermal buckling 
loads. Therefore, equation (22) needs to be rewritten in terms 
of temperature rather than load factors. For the uniform tem- 
perature case, the thermal forces have the following forms: 

N x = (.Ana* + A 12 a y ) T 
N y = (i42lQ* + i422<*y) T 
Nj y = A 66 a xy T (28) 


which were obtained from equations (4) and ( 6 ) setting T\ = 
T 2 = T. 


Case 3. 2C2S edge condition 
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The coefficient o“ nJW appearing in equation (23) contains 
thermal forcing terms (ref. 29). Thus, may be written 

in two parts as 


+ [e(m, k)Nj + C (n, *)<] (29) 


where is the first part of a]^ nki without the thermal 

forcing terms, £(m, k) and C(n,^) are » respectively, the nu- 
merical coefficients of Nj and N y , whose values change with 
the indicial and the edge conditions. 


In light of equation (29), equation (22) could be rewritten 


as 
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where Af mn « is the modified M mn kt in which o” nW is re- 
placed with Timnkti and 


+ C(n,e)(A 2 ia x + A 22 a v ) 


(31) 


In equation (30), both Mmnkt P mnkt terms contain 

material properties that are temperature dependent. Thus, 
in the eigenvalue solution process using equation (30), one 
has to assume a temperature T a and use the correspond- 
ing material properties as inputs to calculate the eigenvalue 
1/Tcr where Ter is the buckling temperature. This iteration 
process must be continued until the assumed temperature T a 
approaches the buckling temperature Ter- 

Thus, in thermal buckling, the eigenvalue solution process 
requires a temperature iteration process and, therefore, is 
slightly different from that in mechanical buckling for which 
only a one-step eigenvalue solution process is needed. 

Different Face-Sheet Temperatures (Nonzero 
Thermal Moments) 

When face-sheet temperatures are different (i.e., T\ ^ T 2 ), 
the sandwich panel will be subjected not only to ther- 
mal forces but also to thermal moments 

AfJ y }. The problem then becomes a bending one 
and no longer an eigenvalue problem. The panel deflection 
w can then be calculated in terms of Fourier coefficient A mn . 
The buckling condition will correspond to that when the term 
in series representation of w (associated with a particular 
buckling mode shape) becomes unbounded (i.e., A mn — 1 ► oo 
for a given {m, n}). For this case the 4S edge condition will 
be analyzed as an example. 


Thermal Moments 


Let the thermal moments {AfJ, be expressed 

in double Fourier series in accordance with the deformation 
functions given in equations (9) to (11) as 
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In light of the deformation functions (eqs. (9) to (11)) and 
the thermal moment expressions (eqs. (32) to (34)), the en- 
ergy equations (7) (setting u = v = 0) and (8) may be inte- 
grated to yield the forms given in reference 29. The thermal 
moment terms turned out to be linear functions of A mn , B mn , 
or C m „, and not quadratic functions of { A mn , B mn , C mn } in 
the energy expressions (ref. 29). 


Nonhomogeneous Equations 

After the application of the Rayleigh-Ritz method accord- 
ing to equation (21), one obtains three nonhomogeneous si- 
multaneous equations for each indicial set of {m,n}, which 
may be combined to yield one equation in terms of Au 
(ref. 29): 
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where M mnmn is defined in equation (23) under the con- 
ditions m=fc, n=l t r/=32, and the nonhomogeneous term 
Rmnmn is defined by 
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„33 


_ r ,,23 \ 


▼ ^mnmn {K mn 

u mnmn 

755 -33 _ -53 735 

u mnmn u mnmn '*mnmn mnmn 
, ^mnmn(^ mn ^mnmn ) 

+ -a? 155 


(39) 


where J mn , K mn , and L mn are 

defined as 

Jmn — P mn ^ 

=)’+*-( 




nir\ 

Kmn = P rnn 

(— )+5 mn ( 

t) 


/mn\ 

fnirv 

Pmn — Kmn ^ J “h ^mn 1 

It) 



(40) 

(41) 

(42) 


Equation (38) forms an infinite number of nonhomogeneous 
simultaneous equations, each of which is associated with a 
set of {m,n} values (or mode shape) for the calculation of an 
infinite number of Fburier (or Ritz) coefficients A mn in the 
series representation of panel deflection w(x, y) (eq. (9)). 


Buckling Condition 


where the Fourier coefficients K mn i Smn are given by 

4 /■“ f b ,/T . mira . rmry , . 

Fmn= riJ J 0 M * sm T" sin ( 35 ) 

4 r f b 1#T . mirx . miry . 

Hmn = ab Jo Jo My sin ~ sin ~r dxdy (36) 


The calculated Ritz coefficients A mn 
functional form: 


Amn — 


[ ]r 


have the following 
(43) 
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where the numerator [ ] mn contains {Af mnmni ftrmmn, 
and the denominator A is the determinant of the 
coefficients of Aiu of the nonhomogeneous simultaneous equa- 
tions written out from equation (38). The detailed expres- 
sions of A of order 12 is shown in references 9 through 11 
and 29 for either symmetrical or antisymmetrical buckling. 

The mathematical meaning of the buckling state in light 
of equation ( 43 ) is that the Ritz coefficient A mn becomes 
unbounded (i.e., infinite panel deflection, or A — ► 0). That 
is, when the budding state is reached, the term in the series 
(eq. ( 9 )) that corresponds to the particular deformation mode 
shape becomes the most important term. 

From the above analysis, one sees that the budding condi- 
tions for the cases with and without the thermal moments are 
identical (i.e., A = 0) under the classical small deformation 
theory. Because of this finding, similar bending analyses for 
nonzero thermal moments for other edge conditions were not 
carried out. 

NUMERICAL EXAMPLES 
Physical Properties of Sandwich Panels 

The sandwich panel is assumed to be fabricated with tita- 
nium face sheets and titanium honeycomb core, having the 
following geometrical and material properties: 

Geometry: h = 1.2 in.; t M = 0.032 in.; a = a 0 = 24 in. (for 
varying 6 ), or o 6 = oj (for constant panel area) 

Material properties: 


Face sheets 



70 °F 

1000 °F 

E x = E y , lb/in 2 

16 x 10 6 

10.5 x 10 6 

G sy , lb/in 2 

6.2 x 10 ® 

4.7 x 10 6 

V xy = Vyx 

0.31 

0.31 

a x = ay, in/in-°F 

4.85 x 10 " 6 

5.6 x 10 - 6 

Ox v , in/in-°F 

0 

0 


Honeycomb core (properties at 600 °F) 


G Cyx = 0.81967 x 10 5 lb/in 2 
Gc*x = 1-81 x 10 5 lb/in 2 


Buckling Interaction Curves 

In generating the data for plotting the budding interaction 
curves, ky was set to zero. For a given different values of 
k x were assigned; then the corresponding eigenvalues l/fc* y 
were calculated from equation (22). Typical buckling interac- 
tion curves plotted in kx -k xy space for square panel (£ = 1 ) 
are shown in figure 4. The additional set of buckling interac- 
tion curves shown in broken curves is for the case when the 
effect of transverse shear is neglected. Fbr the square panels, 
the buckling interaction curves for 4S, 2S2C, and 4C cases 
are continuous curves of symmetric buckling. However, for 
the 2C2S case, the buckling interaction curves are composite 


curves, partly for symmetric buckling and partly for anti- 
symmetric buckling. Notice that the effect of the transverse 
shear is quite large for the sandwich panel. Without a con- 
sideration of the transverse shear effect, the sandwich panel 
buckling strength could be overpredicted considerably. The 
4S case has the lowest buckling strength. Through clamp- 
ing two opposite edges (i.e., from 4S case to 2 C 2 S and 2S2C 
cases), the buckling strength could be enhanced considerably. 
By additional damping of the other two opposite edges (i.e., 
from 2C2S and 2S2C cases to 4C case), further improvement 
of the buckling strength could be achieved. However, the im- 
provement is not as large as that for the previous case (i.e., 
from 4 S case to 2C2S and 2 S 2 C cases). With or without a 
consideration of the transverse shear effect, the improvement 
of buckling strength through edge dampings is larger in pure 
compression than in pure shear. 



Fig. 4. Combined load budding interaction plots for a 
honeycomb-core sandwich panel under different edge condi- 
tions. 


Buckling Interaction Surfaces 

By using families of buckling interaction curves gener- 
ated for different values of ^ (a = constant) and for differ- 
ent edge conditions, three-dimensional buckling surfaces were 
constructed in [fc*, fc* v , ||] space as shown in figure 5. In the 
figure, the domains of symmetric and antisymmetric buck- 
ling (lowest buckling modes) are also shown. Figure 5 shows 
better visualization of the buckling behavior of the sandwich 
panel than the traditional buckling plots of fc* as a function 
of and k X y as a function of For slender rectangular 

panels (i.e., ^ < 1 ), antisymmetric buddings occur mostly in 
the compression-dominated regions. For wider panels (i.e., 
£ > 1 ), the antisymmetric buddings take place in the shear- 

dominated regions. In the neighborhood of 3 = 1, the lowest 
buckling modes are all symmetric (i.e., m = 1 , n = 1 ) for 
the 4S, 4 C and 2 S 2 C cases, and only for the 2C2S case, the 
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lowest buckling mode in the compression- dominated region is 
antisymmetric (i.e., m = 2, n = 1). Such buckling behavior 

also occurs in the flat rectangular plates of ^ « 1 . 




Fig. 5. Buckling interaction surfaces for honeycomb-core 
sandwich panels under different edge conditions (a = con- 
stant). 



(d) Two sides simply supported, two ends clamped (2S2C). 
Fig. 5. Concluded. 


Buckling Curves for Pure Compression and Pure 
Shear 

The buckling surfaces shown in figure 5 were constructed 
under the condition a = constant, and may not serve as ideal 
design plots for aerospace structural panels, because, when 
h. js changed (under a = constant), the panel weight (i.e., 
panel area afc) also is changed accordingly. In the aerospace 
structural designs, the main objective is the structural op- 
timization. That is, for a given panel weight, the objective 
is to search for a panel with optimum buckling strengths (or 
stiffnesses). For this reason, modified buckling load factors 
k x and fc*y (ifcv = 0 ) were recalculated as functions of 5 
under the condition 06 = 0* = constant (instead of o = con- 
stant). Figures 6 and 7, respectively, show the alternative 
plots of kx as a function of ^ for pure compression and k xy 

as a function of 3 for pure shear when the panel area ah 
was kept unchanged. In practical applications, the panel has 
to be supported by an edge frame (cross section assumed 
constant); therefore, the edge frame weight (or edge frame 
length, (o + b)/2a 0 ) was also plotted in figures 6 and 7 as 
a function of The square panel (5 = 1) has the mini- 
mum edge frame weight; however, it does not have the opti- 
mum buckling strengths both in compression and shear. The 
compressive budding strengths (fig. 6 ) reached minimum at 
£ = 1 . 6 , 1 . 4 , 2.2, and 1 . 0 , respectively, for the 4S, 4C, 2 C 2 S, 
and 2S2C cases. The lowest shear buckling strengths (fig. 
7) occur at 5 = 0.9, 0.9, 1.2, and 0.7, respectively, for the 
4 S, 4C, 2C2S, and 2S2C cases. Figures 6 and 7 serve as de- 
sign curves for selecting the desired sandwich panel geometry 
(i.e., | value). To boost the panel budding strengths both in 
compression and shear, some weight penalty resulting from 
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the edge frame is inevitable. The desirable high stiffness-to- 
weight-ratio panel shapes will be slightly slender < 1), 
especially in light of the compressive buckling strength. 



Fig. 6. Comparison of compressive buckling strengths of 
honeycomb-core sandwich panels under different edge con- 
ditions (constant panel areas). 



Fig. 7. Comparison of shear buckling strengths of 
honeycomb-core sandwich panels under different edge con- 
ditions (constant panel areas). 

Thermal Buckling Curves 

For most of the practical materials, the coupling coefficient 
of thermal expansion a xy is zero. Therefore, in generating the 
data for thermal buckling curves, the following conditions 
were imposed: a x — a y , = 0 (i.e., Nl y — Mj y = 

0), Ml = Mj - 0. 

These conditions will induce in-plane biaxial thermal com- 
pression without shear and bending. Figure 8 shows the 
buckling temperature, 7^, plotted as a function of with 
the panel length o kept constant. Those thermal (biaxial 


compression) buckling curves somewhat resemble the uniax- 
ial compressive buckling (mechanical buckling) curves k x as 
a function of ^ shown in figure 5. For buckling temperatures 
higher than 1000 °F, the face sheet material property data at 
1000 °F were used as inputs to equation (30) for calcu- 
lations because of the lack of material property data at high 
temperatures. For the honeycomb core, the only available 
material property data at 600 °F had to be used as inputs 
for Ter calculations. The buckling temperature was found 
to be relatively insensitive to the material property change 
with temperature. 

For the present particular panel (i.e., dimensions chosen), 
the thermal buckling temperatures Ter exceed the titanium 
melting point (3074 °F) at the low ^ and gradually decrease 

with the increase of At high aspect ratios, Ter for 4S 
and 2C2S cases level off at about 1000 °F (below superplas- 
tic temperature, 1650 °F), and for 4C and 2S2C cases, at 
temperatures slightly below the melting point. 


9 x 10 3 



Fig. 8. Thermal buckling temperatures for honeycomb-core 
sandwich panels under different edge conditions (a = con- 
stant). 

Figure 9 shows the alternative plots of as a function 
of ^ for constant-area panels (i.e., ab — constant). The 
lowest buckling temperatures for 4S, 4C, 2C2S, and 2S2C 
cases are, respectively, at 1297 °F, 3702 °F (above melt- 
ing point), 2194 °F (above superplastic temperature), and 
2205 °F (above superplastic temperature) and occur, respec- 
tively, at — = 1.0, 0.975, 1.8, and 0.5. 

For the present sandwich panel, the actual thermal buck- 
ling will take place only for the 4S case in the region 1.5 
< ^ < 1.8. Outside this region for the 4S case and for all 
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the range of £ for the other three edge conditions, no actual 
thermal buckling could occur because the sandwich panel will 
first undergo superplastic creep or melting. 



Fig. 9. Thermal buckling temperatures for honeycomb-core 
sandwich panels under different edge conditions (constant 
panel areas). 


CONCLUDING REMARKS 

By using the Rayleigh-Ritz method of minimizing the to- 
tal potential energy of a structural system, the combined load 
(mechanical or thermal) buckling equations were established 
for orthotropic rectangular sandwich panels supported under 
four different edge conditions. Two-dimensional buckling in- 
teraction curves and three-dimensional buckling interaction 
surfaces were constructed for high- temperature honeycomb- 
core sandwich panels. The buckling interaction surfaces pro- 
vide easy visualization of the variation of the panel buckling 
strengths and the domains of buckling modes (symmetric and 
antisymmetric) with the edge condition. Furthermore, the 
buckling temperature curves for the sandwich panels were 
presented. 

The effect of transverse shear on the buckling strength is 
quite large for sandwich panels, and by neglecting the trans- 
verse shear effect, the buckling strengths could be overpre- 
dicted considerably. With the inclusion of the transverse 
shear effect, the buckling load factors became panel-size de- 
pendent in addition to panel-aspect-ratio dependent. Clamp- 
ing the edges could enhance the buckling strength greatly 
more in compression than in shear. Thermal buckling condi- 
tions for the cases with and without thermal moments were 
found to be identical for the small deformation theory.* 

*The author gratefully acknowledges the contributions by Barry 
Randall in setting up computer programs for the eigenvalue 
extractions. 
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